Multiscale nature of hysteretic phenomena: Application to CoPt-type magnets 
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We suggest a workable approach for the description of multiscale magnetization reversal phe- 
nomena in nanoscale magnets and apply it to CoPt-type alloys. We show that their hysteretic 
properties are governed by two effects originating at different length scales: a peculiar splitting of 
domain walls and their strong pinning at antiphase boundaries. We emphasize that such multiscale 
nature of hysteretic phenomena is a generic feature of nanoscale magnetic materials. 



The problem of magnetization reversal (MR) is a cen- 
tral and longstanding problem in the theory of mag- 
netism. Main complications in its treatment stem from 
the fact that effects due to magnetostatic and elastic 
forces manifest themselves in the whole nanoscale range 
while local interactions must still be treated on the atom- 
istic scale. As the length scale increases, new physical 
effects arise, but coarse-graining over the lower-scale de- 
grees of freedom is not easily accomplished. Specifically, 
MR is realized by motion of domain walls (DW), the two- 
dimensional surfaces where magnetization changes its di- 
rection. DW's may easily move unless they are pinned by 
defects, which implies that MR properties are not intrin- 
sic and depend strongly on the microstructure. Exper- 
imentally, this results in huge variations of the coercive 
force H c and other properties depending on material pro- 
cessing. It is crucial to address the DW structure and its 
interaction with microstructure on different length scales 
ranging from interatomic to submicron. 

In modern micromagnetic methods jjj complex mi- 
crostructures of hard magnets are defined by (usually 
piecewise constant) fields of macroscopic material prop- 
erties such as the exchange constant A(r) and the mag- 
netocrystalline anisotropy (MCA) constant K(r). Such 
calculations were mainly focussed on the role of grain 
boundaries; since they are extremely hard to describe mi- 
croscopically, their properties are postulated in an ad hoc 
manner making reliable predictions problematic. From 
the other hand, typical microstructures of hard mag- 
nets with Llo crystal structure including CoPt, FePt and 
FePd are dominated by other defects g §, [§ §, §, §, 
antiphase boundaries (APB) and twin boundaries (TB), 
which are crystallographically coherent and easier to 
treat than grain boundaries. Although these materials 
were used and studied for many decades and the forma- 
tion of their microstructure is well understood, the mech- 
anism of MR is still a mystery. The goal of the present 
letter is to lay the foundation of a consistent multiscale 
theory of MR in CoPt type magnets. 

Microstructural evolution during Llo ordering is 
strongly affected by tetragonal lattice distortions J|, |). 
After a relatively short 'tweed' stage of annealing after 
quench, when the ordered domains achieve some char- 
acteristic size lo ~ 10 nm, they develop so-called poly- 
twinning, i.e. the formation of regular arrays ('stacks') 



of ordered bands ('c-domains') separated by TB's 
|, U | §. The tetragonal axes c of the c-domains 
(pointing along one of the three cubic axes of the par- 
ent fee phase) alternate regularly making 90° angles be- 
tween the adjacent domains (below we use the obvi- 
ous terms 'X-domain', 'XY-stack', etc.). In addition to 
polytwinning, the observed and simulated microstruc- 
tures always contain a high density of APB's in the c- 
domains 1,111,0. 

Polytwinned magnets were extensively studied with 
analytical micromagnetic methods §. Due to high MCA 
the DW width 5 = tt^A/K is quite small (5-10 nm), 
while the anisotropy field H a = 2K/M (M is the sat- 
uration magnetization) significantly exceeds the typical 
magnetostatic field H m = 4ttM (parameter r\ — H m /H a 
is close to 0.1 in CoPt and FePt and 0.38 in FePd §). 
Therefore, at d > 5 (d is the c-domain thickness) each 
c-domain may be regarded as an individual magnetic do- 
main with intrinsic 90° DW's at the TB's g. It is as- 
sumed that MR is associated with 'macrodomain walls' 
(MDW) crossing many c-domains. Such MDW's were ob- 
served experimentally g, ^] , but their internal structure 
and mechanisms of coercivity are unknown 0]. 

Microscopic mean-field method. Consider a binary al- 
loy AB with the classical Hamiltonian 



H = H con{ {n,} + ^2 n i n j \-JijfiiPj + A'A. 

i<j 



■22 n i[ e i(Pi) ~ HoA*i] (1) 



where -ff C onf is the configurational part of the Hamilto- 
nian; i and j run over lattice sites; m = 1 if site i is 
occupied by a magnetic atom A and rii = otherwise 
(let B be non-magnetic for simplicity); pi is the rigid 
classical magnetic moment of the atom at site i; Jij, the 
Heisenberg exchange parameters; Ho, the external mag- 
netic field; ei(fii), the MCA energy equal to —biOt^i) 2 
for easy-axis anisotropy; and Dij , the magnetic dipole- 
dipole interaction tensor. 

The free energy in the mean-field approximation 
(MFA) may be related to the 'mean fields' H; = Ho + 



- Dij)cjm.j where m t = (/Zj) and 



are 



2 



local magnetizations and concentrations, respectively: 



F : 



-Ej^ DD -T^Cjbx J d% exp [^(H^ - e»)] . (2) 



Dij)m.j is the total ex- 



Here E JiDD — J2 c i c ] m i(-Ji] 

change and dipole-dipole energy. Equilibrium states may 
be found using self-consistency relation nij = —dF/dtli. 

If magnetization M(r) slowly varies in space and is 
constant in magnitude, Eq. (Q) reduces to the micromag- 
netic free energy ||. In this case all choices of Jy and 
6j in the defect-free regions are equivalent if they pro- 
duce the same macroscopic properties A and K. How- 
ever, variation of Jy and 6$ near defects like APB's must 
be studied using first-principles techniques. Microscopic 
MFA calculations with these parameters may be used to 
describe DW interaction with a defect at the length scale 
of S. At larger, microstructural length scales micromag- 
netic methods jj| may be used with singularities of A 
and K at the defects. However, in hard magnets the mi- 
croscopic approach also turns out to be convenient for 
the studies of regions containing up to ~ 10 6 atoms; in 
such calculations some model and bi reproducing the 
actual defect properties may be used. 

Below, we first address the equilibrium structure of 
MDW's in a defect-free system and find a peculiar MDW 
splitting phenomenon. Then we explore the behavior of 
MDW's in external field, briefly report our results for 
pinning at APB's, and finally deduce the combined effect 
of MDW splitting and pinning on coercivity. 

Equilibrium macrodomain walls. It is convenient to 
use the body-centered tetragonal (bet) representation of 
the fee lattice (with c/a = y/2 and c equal to the fee 
lattice parameter). Two opposite edges of a rectangular 
simulation box are aligned with the (110) TB's in the 
XY-stack; MDW's are normal to them. The boundary 
conditions are periodic, and the dipole-dipole fields are 
computed using the Fourier transforms. 

Equilibrium MDW's of two characteristic (110) and 
(001) orientations are shown in Fig. [l] for the CoPt model. 
Infinite stack of c-domains is assumed with ideal TB's in 
the (110) planes. Fully ordered c-domains have the same 
thickness d = 64a ~ 17 nm, which corresponds to an 
early stage of annealing shortly after the polytwinned 
stacks are formed |6) . The anisotropy is uniaxial with 

pointing along the local direction of the tetragonal 
axis and bi = b. Room temperature T = 0.4T C is as- 
sumed averywhere (for CoPt T c ~ 720 K). For sim- 
plicity, only 3d-metal atoms are assumed to have mag- 
netic moment fi. Parameters b and fi were chosen so 
that MFA gives experimental room-temperature values 
of K = 4.9 x 10 7 erg/cm 3 and jj = 0.082 for CoPt @. 
The parameters Jy for nearest and next-nearest neigh- 
bors were chosen as J2/J1 = 2/3, J3/J1 = 1/6; A is 
isotropic with this choice. 

MDW's of both orientations shown in Fig. [I] have a 
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FIG. 1: Macrodomain walls in an ideal XY-stack with 
d = 64a bc t oriented normal to: (a),(b) (HO); (c),(d) (001). 
Frames (a),(c): magnetization M(r) (arrows and color); 
(b),(d): dipole fields (arrows) and magnetic charge density 
p = — div M (color). The axes at top refer to (a),(b); in (c),(d) 
they are rotated 90° around the horizontal line. Sticks with 
circles show vector directions for cells i = (8ni +4, 8ri2 +4, 1) 
with integer m, ri2- The length of a stick (up to the center 
of the circle) is proportional to the vector projection onto the 
graph plane; the diameter of full (open) circles, to the positive 
(negative) out-of-plane component. Small points at the end of 
long sticks show their direction. Color in (a),(c) (guide to the 
eye) is obtained by mixing red, green and blue according to 
the values of three vector components. In (b),(d) the intensity 
of red (blue) color is proportional to the relative density of 
positive (negative) magnetic charges p/|p ma x|. Twin bound- 
aries are vertical. The simulation boxes had 512x128x1 and 
512x128x2 bet cells for (HO) and (001) MDW's, respectively. 
Figures are 2 times higher than the simulation boxes; they 
are trimmed at top and bottom to conserve space. 



peculiar feature which turns out to be the key to the 
understanding of MR in polytwinned magnets: the DW 
segments (DWS) located in adjacent c-domains are dis- 
placed in respect to each other. Such configuration 
with regularly alternating displacements is beneficial for 
the exchange energy because magnetization within each 
DWS is parallel to that in the adjacent c-domains (note 
that the color of horizontal DWS's in Fig. |l|a,c matches 
that of adjacent c-domains). The magnetostatic energy 
for this MDW configuration is higher compared to the 
rectilinear one due to the presence of short magneti- 
cally charged (in micromagnetic terms Q) segments of 
90° DW's at the TB's (note the red and blue 'blobs' in 
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Fig. [j]d). Actual MDW configuration emerges as a re- 
sult of competition between these interactions. However, 
due to small rj the dipole-dipole interaction is, in fact, 
unimportant at the length scale of 5. 

As it is clear from Fig. [l], the magnetic charges both 
at the DWS's of a (lTO) MDW and at the TB's alternate 
in sign, so that each MDW is magnetically uncharged as 
a whole, and its magnetic field quickly falls off at r > d. 
DWS's of (001) MDW's do not carry magnetic charge. 

The orientation of DWS's within c-domains in (110) 
MDW is determined by the anisotropy of the exchange 
constant a a = A±/Au (A± and An are values of A nor- 
mal and parallel to c), and by the parameter £ = rjd/8. 
At £ <C 1 the DW orientation is determined solely by cxa- 
At £ 3> 1 the DW aligns parallel to the tetragonal axis, 
and ctA is irrelevant. The (110) MDW shown in Fig. [j]a,b 
is in the crossover region with £ ~ 0.2 and a a — 1. 

Partial macrodomain walls. MDW's behave remark- 
ably in external magnetic field Ho because DWS's are 
held together only by magnetostatic forces. If these forces 
were absent (77 — ^ limit), each DWS would be able to 
move freely until it meets another DWS in an adjacent 
c-domain. Some of this freedom remains at small 77. Con- 
sider an XY-stack with (110) TB's , as above. Compo- 
nent H z does not exert any force on DWS's, and the effect 
of Ho on MDW's depends only on the orientation of its 
projection Hi on the xOy plane. If H_l is normal to the 
TB's (H x — H y ), then the forces acting on all DWS's are 
in the same direction, and both MDW's shown in Fig. |l| 
move as a whole like usual DW's in a homogeneous crys- 
tal. However, if H^ is parallel to the TB's (H x — —H y ), 
the forces acting on X- and Y-DWS's differ in sign, and 
the total force acting on the MDW is zero. 

Let us call a set of DWS's of a MDW in all even or 
odd c-domains as a partial macrodomain wall (PMDW), 
e.g. X-PMDW. There are two MDW realizations dif- 
fering in the sign of the relative displacement L of the 
two PMDW's. For one of them Hj_ increases \L\ ('split- 
ting'), and for the other, reduces it ('swapping'). If H± is 
inverted, splitting turns to swapping and vice-versa. Ac- 
cordingly, two threshold fields may be introduced. Let 
us fix, for instance, Y-PMDW and examine the effect of 
H x on X-PMDW. The threshold splitting field H sp is the 
minimal value of H x required to move the X-PMDW far 
from the Y-PMDW. The threshold swapping field H SVJ is 
required to 'drag' the X-PMDW through Y-PMDW in 
the opposite direction, overcoming the exchange barriers 
at the TB's . H sw is of lesser importance than H sp be- 
cause even if _ff sp < H sw , the existence of some MDW's 
in the 'splitting relation' with H x is sufficient to com- 
plete the MR. Note also that H sw cx <i _1 and quickly 
falls off with increasing d. Both H sp and _ff sw are 'intrin- 
sic' properties for a given d, i.e. they exist even if there 
are no defects (APB's) in the c-domains. 

H sp may be easily estimated micromagnetically for a 
(001) MDW. When its two PMDW's are displaced at dis- 
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FIG. 2: Split (001) macrodomain wall in a stack with d 
32abct- Presentation is similar to Fig. Id. 



tance L, the alternating magnetic charges appearing at 
the TB's form a regular stack of 'capacitors' (see Fig. ||). 
For L > d, the magnetic field is concentrated inside these 
capacitors, and the total magnetostatic energy is propor- 
tional to L; the coefficient in this proportionality deter- 
mines H sp . Such calculation gives H sp = nM (or, al- 
ternatively, H sp = H a r]/4). The splitting threshold for 
(110) MDW is larger because of the magnetic charges on 
its segments. In the above estimate we neglected the de- 
viation of M from the easy axis at distances ~ S from 
the TB's , which reduces H sp for small d. 

Thus, at 77 -C 1 we have H sp <C H a . In particular, 
using the room-temperature values of M |2| , we obtain 
H sp of 2.5 kOc for CoPt and 3.5 kOe for FcPt and FePd. 
These values are close to the typical fields at which hys- 
teresis is observed in these magnets. Quite naturally, 
they are much less than the coercive fields corresponding 
to uniform rotation of magnetization in the c-domains [0 . 
A peculiar feature of a PMDW in the polytwinned crystal 
is that it separates qualitatively different magnetic do- 
mains with different free energies, and hence there is an 
intrinsic force acting on an isolated PMDW (its density 
per unit cross-section of the c-domain is f s — 2MH sp ). 

Domain wall segments, pinning and coercivity. Real 
polytwinned alloys, as we noted above, usually contain 
a high density of APB's in the c-domains. It is known 
that DW's may be pinned by APB's jl0| ], and it was 
argued that this mechanism may explain high coercivity 
of CoPt type alloys 0, || . Here we show that this pinning 
is indeed quite strong, but the MDW 'dissociation' into 
segments also has a profound effect on MR. 

The APB pattern within c-domains determines the sta- 
tistical properties of the random value U(x) where x is 
the coordinate of a DWS, and U is its APB-induced ex- 
cess free energy. For the following basic consideration 
we assume that U(x) describes a distribution of simi- 
lar pinning centers with the typical distance Id between 
them. The maximal slope of U(x) determines the un- 
pinning threshold H u , i.e. the value of the field Ho || c 
required to unpin a single c-DWS. Thus, in the absence of 
magnetostatic forces at H < H u each c-DWS is pinned, 
and at Hq > H u it can move freely to the surface of 
the polytwinned stack. These assumptions are obvious 
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for patterns with isolated APB's as in Fig. 5 of Ref. || 
where Id ~ 5-10d; however, main features of more com- 
plex patterns [Q, M, [| || may also be described by some 
characteristic values of H u and Id ^ d. 

The maximum possible H u is achieved when DWS's are 
parallel to isolated plain APB's. To evaluate this maxi- 
mum in the studied alloys, we explored the modification 
of exchange and anisotropy at an isolated (lOl)-oriented 
APB in CoPt, FePt and FePd using the tight-binding 
linear muffin-tin orbital method. We found that MCA is 
strongly suppressed at the APB leading to a DW attrac- 
tion to the APB with H u of about 11,7 and 1.5 kOe for 
CoPt, FePt and FePd, respectively. These values signif- 
icantly exceed the observed coercivities ||. Since high 
MCA is associated with Llo ordering, it is natural that 
local disorder at an APB suppresses MCA. 

Let us now explore the role of intrinsic magnetostatic 
forces in MR. We restrict ourselves to the case of H || Ox 
in a single crystal with different types of polytwinned 
stacks ('single crystal' refers to the parent fee lattice). 
Such field does not affect YZ-stacks, except for a re- 
versible transverse magnetization. In other stacks Ho 
exerts a force density / = 2MHq on X-DWS's only. 

From our definition of H sp it follows that if c-PMDW 
in Xc-stack (c =Y or Z) is held in place (e.g., by pinning), 
then the free (unpinned) X-PMDW moves to infinity at 
H x > H sp (neglecting the demagnetizing effects). At 
H x < H sp the X-PMDW moves to a finite distance L(H X ) 
from the c-PMDW. By definition, L(H sp ) — oo. Except 
for a close vicinity of _ff sp , L < d. 

The path of MR differs qualitatively in two cases: (a) 
H sp > H u , and (b) H sp < H u . In case (a) the MDW 
can only move as a whole (intrinsic forces f s are strong 
enough to unpin a PMDW), but DWS's may to some ex- 
tent 'adapt' to the local pinning potential. Let us define 
the characteristic displacement l m 'allowed' by magneto- 
static forces as l m = L(H U ). Obviously, l m increases with 
H u , and l m — > oo at H u — > H sp . If l m > Id, the effective 
unpinning field of the MDW is Hjj ~ 2H U (Ho affects 
one half of the DWS's, while all DWS's can be effectively 
pinned). However, as l m becomes smaller than Id, Hjj 
is quickly reduced, because only a fraction of DWS's can 
be pinned simultaneously (MDW is effectively 'rigid' and 
can not make 'kinks' comparable with Id). 

Now, in case (b) the MDW may split in two PMDW's 
at sufficiently large Ho. The MDW is 'soft', and each 
DWS can be effectively pinned. During the whole cy- 
cle of remagnetization pinned Y- and Z-PMDW's do not 
move at all, and hence the intrinsic (magnetostatic) force 
acting on a moving X-PMDW will change its sign every 
time that it crosses, e.g., an Y-PMDW (in XY-stack). 
Thus, for an X-PMDW we have H v = H u + H sp . Dif- 
ferent MR properties may, in principle, be observed if 
the sample was previously magnetized to saturation by 
Ho along, e.g., (111). In this case there are no Y- and 
Z-PMDW's, and in some stacks the switching threshold 



is Hjj = H u — H sp . The hysteresis loop in this case may 
have a peculiar double-step feature with magnetization 
jumps at Hq = H u ± H sp . Such 'hysteresis memory' ef- 
fects provide a test of the H sp < H u relation and may be 
used to design 'programmable' magnets. 

Thus, for an initially demagnetized sample with the 
given H u the coercive force H c first rises with H sp , 
reaches its maximum ~ 2H U at H sp ss H u , stays roughly 
constant in the range oo > l m > Id, and finally falls off at 
l m < Id- High coercivity of CoPt is probably due to the 
fact that this magnet with suitable processing is close to 
the optimal H sp ~ H u condition. By contrast, in FePd 
H sp is too large (H sp / H u w 2.3 even for maximal possi- 
ble H u ), so that l m <§; d; this conclusion agrees with the 
relatively small observed H c /H a ratio. 

In our analysis we neglected the demagnetizing fields. 
Close to saturation they are comparable to H sp , but they 
have little effect on H c , because they are small when the 
total magnetization is close to zero. 

In conclusion, we developed a workable technique for 
the description of multiscale MR phenomena in hard 
magnets and applied it to CoPt-type alloys. We have 
shown that coercivity of these materials has at least two 
sources originating at different length scales: strong pin- 
ning of MDWs by APB's and their splitting at TB's . 
This DW splitting seems to be a generic effect that may 
also affect MR in other groups of materials, including 
some nanostripes |ll| and multilayers; it may also occur 
at grain boundaries of certain misorientations. 
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